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Abstract 

Hybrid adiabatic potentials are considered in the framework of the QCD string 
model. The einbein field formalism is applied to obtain the large-distance behaviour 
of adiabatic potentials. The calculated excitation curves are shown to be the result 
of interplay between potential-type longitudinal and string-type transverse vibrations. 
The results are compared with recent lattice data. 

1 Introduction 

There is no doubts now that hadrons with explicit gluonic degrees of freedom should ex- 
ist. This idea is supported not only by general arguments from QCD, but also by lattice 
simulations of pure Yang-Mills theory. In the absence of exact analytical methods of nonper- 
turbative QCD one relies upon models to describe gluonic mesons, so that the challenging 
question arises of how to introduce effective degrees of freedom for soft (constituent) glue. 

There is a lot of indications now that gluonic mesons are already found experimentally, 
but the conclusive evidences have never been presented; there is no hope that in the nearest 
future data analyses could shed light on this problem and to offer necessary feedback for 
model building. The current situation is such that the predictions of different models on 
hadronic spectra and decays are involved in order to pin-point the signatures for gluonic 
mesons. 

On the other hand, lattice calculations are now accurate enough to provide reliable data 
on constituent glue and to check model predictions. In this regard recent measurements of 
gluelump and hybrid adiabatic potentials are of particular interest. These simulations 
measure the spectrum of the glue in the presence of infinitely heavy adjoint source (gluelump) 
and in the presence of static quark and antiquark separated by some distance R. These 
systems are the simplest ones and play the role of hydrogen atom of soft glue studies, as the 
complicated problem of the centre-of-mass motion separation is not relevant here. 

Hybrid adiabatic potentials enter heavy hybrid mass estimations in the Born-Oppenheimer 
approximation: these potentials are to be inserted into QQ Schroedinger equation in order 
to obtain spectra of hybrids with heavy quarks. The large R limit is interesting per se, as 
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the formation of confining string is expected at large distances, so the direct measurements 
of the string fluctuations become available and the possibility exists to discriminate between 
different models of the effective string degrees of freedom. 

2 Constituent gluons at the end of the string 

The notion of confinement is usually described in terms of area law asymptotics for the 
Wilson loop expectation value, defined as an integral along some closed contour C, averaged 
over gluonic vacuum configurations: 

(W(C)) = Tr(Pexp ig <f A,dz,), (1) 

where trace is taken over colour indices. The area law asymptotics implies that 

(W(C))^N c exp(-aS), (2) 

where Nc is the number of colours, a is the string tension, and S is the surface bound by 
the closed contour C. As the initial expression ([]]) depends only on the contour, the area in 
(0) should depend only on the contour too, and should be the minimal area. 

The area law asymptotics provides the action of the string, and in the case of minimal 
area this string is also "minimal". The effective string model should be arranged to allow the 
extra degrees of freedom to populate the string and to be responsible for more complicated 
string configurations. In what follows these extra degrees of freedom are defined in the 
framework of the QCD string model. This model deals with quarks and point-like gluons 
propagating in the confining QCD vacuum, and is based on Vacuum Background Correlators 
method 0. 

The QCD string model for gluons is derived from the perturbation theory in the nonper- 
turbative confining background developed in The main idea is to split the gauge field 
as 

A tl = B tl + a M , (3) 

which allows to distinguish clearly between confining gluonic field configurations and 
confined valence gluons a M . Confining QCD vacuum is given by the set of gauge invariant 
field strength correlators made of B^, which are responsible for the area law asymptotics 
while the valence gluons are treated as perturbation at this confining background. 

The starting point is the Green function for the gluon propagating in the given external 
field ||: 

G^x, y) = (D 2 (B)5, U + 2igF lxv {B))-\ (4) 

The term, proportional to F^B) is responsible for the gluon spin interaction, it can be 
treated as perturbation ||, and we neglect it for a moment. The next step is to use Feynman- 
Schwinger representation for the quark-antiquark-gluon Green function, which, for the case 
of static quark and antiquark, is reduced to the form 

G(x g ,y g ) = JdsJ Dz g exp(-K g )(W} B , (5) 
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Figure 1: Hybrid Wilson loop. 



where angular brackets mean averaging over background field. The quantity K g is the kinetic 
energy of gluon (to be specified below), and all the dependence on the vacuum background 
field is contained in the generalized Wilson loop W, depicted in Fig.l. 

The main assumption of the QCD string model is the minimal area law, which yields for 
the configuration W the form |5j 

N 2 - 1 

(W) B = exp(-a(S 1 + S 2 )), (6) 

where Si and S 2 are the minimal areas inside the contours formed by quark and gluon and 
antiquark and gluon trajectories correspondingly. 



3 Einbein field form of the gluonic Lagrangian 

To define the form of gluon kinetic energy we note that the action of a particle in the external 
vector field is invariant under reparametrization transformations, and, of course, it remains 
true after averaging over the background. So, to proceed further we are to fix the gauge, 
and the most natural way to do this is to identify the proper time r of the path integral 
representation with the physical time x° g . Then the action of the system can be immediately 
read out of the representation (||): 

-f + T" " * Jo d^iw^y-wfw' 2 - 

-a jT 1 d(3 2 ^{w 2 w' 2 y-wiw'i} , (7) 

where the minimal surfaces S\ and S 2 are parametrized by the coordinates w^r, /3j), i = 
1,2, w ifl = -Tf^, w' ifM = -j^f-- In what follows the straight-line ansatz is chosen for the 
minimal surface: 

R 

w i0 = r, w 1:2 = ±(1 "/?)- + (3r. (8) 

The quantity /i = /i(t) in the expression (|7|) is the so-called einbein field 0; here one is forced 
to introduce it, as it is the only way to obtain the meaningful dynamics for the massless 
particle. 
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Let us introduce another set of einbein fields, z/, = z/j(r, to get rid of Nambu-Goto 
square roots in (|7p |J. The resulting Lagrangian takes the form 



2v x 

JO IV2 Jo 2 



2? 2 2 ), 



l l,2 = r 2 - j-(ri, 2 r) 2 , f lj2 = f± — , (9) 
'1,2 z 

and the corresponding Hamiltonian reads 

h = * + 1 + [ + r «^ + r < + r «at- « io) 

" 2 ■ 



20+ j x + j 2 



2A(/i + Ji + J 2 ) \ rf 
2<7i<7 2 



/(^ Jl(A+Jl) + (^)! j2(/1+j2)+ 



2 2 vrif 2 )(pri)(pr 2 ) [ , (11) 

r l r 2 J 



A = (// + Ji)0 + J 2 )- JiJ 2 ^^, Ji= fw-^A), z = l,2. 



rfr 2 jo 

At first glance, the Hamiltonian ( [TT]) looks tractable. Clearly, quantities /i and Ui play the 
role of gluon constituent mass and energy density distributions along the string respectively. 
Nevertheless, introducing einbeins does not do miracles for us. These redundant variables 
are to be found from the conditions 

dH 5H 

-£- = 0, -T7^T = 0, 12 

as the equations (|T2"D play the role of second class constraints. One should do it before 
quantization and substitute the resulting values into the Hamiltonian. Such procedure is 
hardly possible analytically even at the classical level, and after quantization these extremal 
values of einbeins would become nonlinear operator functions of coordinates and momenta 
with inevitable ordering problems arising. 

In what follows we use the approximate einbein field method, which treats the einbeins 
as c-number variational parameters. The eigenvalues of the Hamiltonian ( [TT]) are found as 
functions of \i and z/j and minimized with respect to eibeins to obtain the physical spectrum. 
Such procedure works surprisingly well in the QCD string model calculations, with the 
accuracy of about 5-10% ||. 
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4 Dynamical regimes of the gluonic Hamiltonian 



Even with the simplifying assumptions of the approximate einbein field method the problem 
remains complicated due to the presence of the terms J, responsible for the string inertia. 
If one neglect these terms, then the einbeins are eliminated explicitly from the Hamiltonian 
(|ITD , and one arrives at the potential model Hamiltonian 

H — Jy? + or\ + 0T2. (13) 
It appears, however, that the neglect of string inertia is justified only for R < 1/a/o" [|10| . 



Indeed, in the einbein field method the potential regime corresponds to the case of Ui inde- 
pendent of Pi. For example, for R <C 1/a/^" one has 

E»m=^{n + 3/2)* + (14) 



/ i„( J R) = 2 1 /V 1 /2( n + 3 /2) 1 /2_ 



aV 2 R 2 



2 5 /2( ra + 3/2)V2' 
(n + 3/2) 1 /V/ 2 3aV 2 R 2 

U h2n{R) ~ 2l/2 + 2 7/2( n + 3/2)V2' 

where n is the number of oscillator quanta. The last line in flTjJ) readily gives Ji t %/ '{i ~ 1/6. 
The situation here is similar to the one in the light quark, glueball and gluelump sectors: 
the corrections due to the string inertia are sizeable, but not large, and can be taken into 
account as perturbation ||. 

The situation changes drastically for the case of large R, R 3> 1/ \fo. Here one has 

4qV 3 (n + 3/2) 2 / 3 aR 

Vn[R) = , Vl,2n{ R ) = ~y- 

In this case = \&R 3> fi n , and the potential regime becomes unadequate. 

The case of large R can be treated exactly in the einbein field method. There are two 
different kinds of excitations, along the QQ axis and in the transverse direction, which are 
decoupled in the limit of large R. In this case one has 

3 CT l/3 ( n 4. 1/9)2/3 9 • 3 1 / 2 
En(R) =*R+ ~ ^t/3 j + -R-^ + A + ^ ( 16 ) 

LR 



where A = is the projection of orbital momentum onto QQ axis (z axis). Note, that 
the sub leading corrections due to the longitudinal and transverse vibrations are different, 
the former behaves as -| 1,/3 , like the pure potential regime (|I5|) , while the latter displays 
string- type behaviour ~ A/R. 

The quasiclassical limit of large A, where only rotations around z axis are taken into 
account, was found for the Hamiltonian (|TTD in fllQfl . The large R limit reads 



E{R) =aR+2V3^-, (17) 
R 
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Figure 2: Adiabatic hybrid potentials in various regimes. Quasiclassical (solid line), potential 
(dashed), and flux-tube (dotted) curves for n z = n p = and A = 1, 2, 3. The lowest curve is 
aR. a = 0.22GeV 2 . 
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Figure 3: Corrections to linear behaviour of potentials. QCD string (solid line), potential 



(dashed), and flux-tube (dotted) curves; n z 



n. 



0; a = 0.22GeV 2 . 



which should be compared with the predictions from naive Nambu-Goto string model 

E(R) = aR + ^- (18) 
R 

in the small oscillation limit. The coefficients tc and 2v3 are close to each other, and 
differ mainly due to the fact that string configurations differ: there are two straight-line 
strings in the QCD string model, and one continuous string in the Nambu-Goto case. The 
quasiclassical excitation curve |I| is shown in Fig.2 together with the Nambu-Goto (|I8~D and 



potential curves. Note the absence of the unphysical divergent 1/R behaviour at small R 
both for quasiclassical and potential curves. 

The dominant subleading regime is defined by longitudinal motion: even if no longitu- 
dinal quanta are excited, there is the contribution of zero longitudinal oscillations in (|16|). 
Still, if the distances are not asymptotically large, the potential regime is substantionally 
contaminated by the string-type transverse vibrations. 

Corresponding corrections to linear behaviour for QCD string (|16|), potential regime (|15|), 
and Nambu-Goto string ( ]T8| ) are shown in Fig. 3. 

There is no QCD string calculations for large distances with proper account of gluon 
spin yet, so the direct comparison with lattice results is premature. Nevertheless, some 
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preliminary conclusions can be drawn. The behaviour (|T6|) displays the most pronounced 
difference between QCD string and other models for constituent glue. In the flux tube 
model JIIJ the string vibrations are caused by string phonons, so one expects the Nambu- 
Goto behaviour fll8D at large distances. In contrast to this, here the string vibrations are 
caused by pointlike valence gluons. In the constituent gluon model with linear potential 
[|l2" | one should have the potential-type behaviour fll5l) , while in the QCD string model the 
confining force follows from the minimal area law, and, as a consequence, the contributions 
from the string inertia leave room for the string-type vibrations. 



5 Full QCD string calculations in the potential regime 

Let us now consider the regime of " small" R relevant to the heavy hybrid mass estimations. 
Actually these distances are not very small: one expects that in the case of very heavy 
quarks the hybrid resides in the bottom of the potential well given by the adiabatic curve. 
The lattice results 0] are not very accurate for small R, but the message is quite clear: the 
bottom of the potential well is somewhere around 0.25 fm for lowest curves. 

If only confining force is taken into account, the QCD string model predicts the oscillator 
potential flT4|) with the minimum at R = 0. However, the minimum is shifted, if the long 
range confining force is augmented by the short range Coulomb interaction, which is taken 
in the form 

V ---^L-l^L + ^L (19) 



The coefficients in (19) are in accordance with the colour content of the QQg system. Note 



that the QQ Coulomb force is repulsive. It is in contrast to the flux tube model ||11|| , where 
the string phonons do not carry colour quantum number, so that the QQ pair is in the colour 
singlet state, and Coulomb interaction is attractive. In the so-called single-bead version of 
the flux tube model the minimum due the confining force is at R = 0, and attractive Coulomb 
force cannot shift it, so the single-bead version seems to be ruled out by lattice data. As 
the gluon energies of [0] lie well below the curves (0), the simple Nambu-Goto regime is 
excluded too. 

Below we present the results of the QCD string model with Coulomb force included, for 
small and intermediate values of R. The calculations were performed in the potential regime, 
and the string inertia was taken into account perturbatively, which is justified by arguments 
given above. 

Varying over /Xj with result = 0X1,2, we obtain the Hamiltonian in the form 

2 2/i K J 2 Vri r 2 J 6R K J 

We calculate unperturbed adiabatic energy levels variationally with Gaussian wave func- 
tions used, 

E jlA {fx,R) = (y jlA (r)\H(^r,R)\V jlA (r)) = (H) jlA , (21) 



= m e c^n* ( ; ) ^ ( 22 ) 
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Figure 4: Hybrid potentials in full QCD string model (thin solid curves) compared to lattice 
results (thick light curves). QQ distance R is in units 2r ~ 1 fm and potentials V are in 
units l/r . 



where s is the spin wave function. The unpertubed adiabatic potentials 

V? lA (R) = E jlA (n*(R),R), 



(23) 
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Table 1: Quantum numbers of levels of Fig. 4 



(a) 


j = l, 1 = 1, A = 0,1 


^ u > n u 


(b) 


3 = 1, 1 = 2, A = 0, 1 




(c) 


3 = 2, 1 = 2, A = 0,1,2 




(d) 


j = 2, 1 = 3, A = 0,1,2 





depend on the extremal value jj* defined from the condition ^Esl^tSl — o. 
The string correction Hamiltonian at v -C ^ after integrating out v reads 

O ( 1 ,o 1 



+ 



# StrinS = -7^(-£? + -^ 
o/i z Vri r 2 



(t r / 1 i 
+ 



6/i 2 IVri ' r 2/ 

2 (r~ ~ r~) { r i PkR i Pk ~ r jPkRkPj + RjP k rjPk ~ RjPkr k Pj]j , (24) 

— * — * 

where L 12 = r 12 xp, L = r x p. 

Here we estimate the contribution of the string correction ([24]) taking it at R = 0, where 
it reads 



V™ = ~ — 9 ' (-);• (25) 
By the same way we consider spin-orbit correction at R = 0, 



V? = -g^>., (26) 

where SX = |(j(j + 1) — /(/ + 1) — 2). We also subtract a constant that corresponds to gluon 
self-energy. From comparison with lattice data 0] its value is found to be £ g i = 473 MeV. 
In Fig. 4 (a)-(d) potential 

VMR) = vU R ) + v i tIins + v $° - S § i (27) 

(solid curves) is compared to lattice potentials from || (thick grey curves). Parameters 
a s = 0.225, a = 0.221 GeV 2 are taken from fit of the lattice Coulomb+linear ground-state 
potential. In the Fig. 4 standard notations from the physics of atomic molecules are used. 

In Table 1 the quantum numbers of levels, shown in Fig. 4 (a)-(d), are listed, in terms of 
j,l, A, and atomic notations. Note, that in the QCD string model the gluon is effectively 
massive and has three polarizations ||, and only two of them are excited with magnetic 
components of field strength correlators, used in the lattice calculations. It is just these 
states which are listed in the Table 1. For more details justifying such correspondance see 
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6 Conclusions 

The first results of the QCD string model for the hybrid adiabatic potentials look rather 
encouraging. We have obtained the reasonable agreement with lattice data at small and 
intermediate interquark distances. The most interesting feature of the QCD string model 
is the large distance behaviour of the adiabatic potentials, with potential-type longitudinal 
and string-type transverse vibrations. The full calculations of large distance behaviour with 
proper account of gluonic spin are in progress now. 
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